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The Effective rfg*g* Vertex at Arbitrary Gluon Virtualities 

A. Ali' [] and A.Ya. Parkhomenko 2 ’ [[ 
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2 Yaroslavl State University, Yaroslavl, Russia 

The effective rf g*g* vertex is studied in the standard (Brodsky-Lepage) and mod¬ 
ified hard scattering approaches for arbitrary gluon virtualities in the time-like and 
space-like regions. The contribution of the gluons in the r/-rrieson is taken into 
account, and the wave-function is constrained using data on the electromagnetic 
transition form factor of the rf meson. Our results have implications for the inclu¬ 
sive decay B —> rf X and exclusive decays, such as B —* rf(K,K*), and in hadronic 
production processes N + N(N) —► rf X. 


I. INTRODUCTION 

The effective coupling involving two gluons and the rf -meson enters in a number of 
production and decay processes. For example, the inclusive decay B —> rfX s [T| and the 
exclusive decay B —> rfK J?. jc|, involve, apart from the matrix elements of the four-quark 
operators, the transitions b —> sg*, followed by g* —> rfg |, i, b —> sgg followed by gg —> 
rf ||> as we ll as the transitions g*g* —> rf and g*g —»■ rf [0. Thus, a reliable determination 
of the vertex function, which is often called the transition form factor, 5(g 2 ,g|,m 2 ,) (here, 
q\ and q\ represent the virtualities of the two gluons) is an essential input in a quantitative 
understanding of these and related decays. Apart from the mentioned 5-decays, the rj'g*g(*> 
vertex plays a role in a large number of processes, among them are the radiative decay 
J/rp — > rf 7 and the hadronic production processes N + N(N) —> rf + A", where N is a 
nucleon. The QCD axial anomaly [|J, responsible for the bulk of the rf mass, normalizes 
the vertex function on the gluon mass-shell, yielding 5(0, 0,m 2 ,). The question that still 
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remains concerns the determination of the vertex for arbitrary time-like and space-like gluon 
virtualities, qf; i — 1,2. A related aspect is to understand the relation between the g'g*g* 
vertex and the wave-function of the //-meson. Stated differently, issues such as the transverse 
momenta of the partons in the 7 /-meson and their impact on the r]'g*g* vertex have to be 
studied quantitatively. 

While information on the g'g*g* vertex is at present both indirect and scarce, its electro¬ 
magnetic counterpart involving the coupling of two photons and the 7 /-meson, P// 7 * 7 , more 
generally the meson-photon transition form factor, has been the subject of intense theoreti¬ 
cal and experimental activity. In particular, the hard scattering approach to transition form 
factors, developed by Brodsky and Lepage 0, has been extensively used in studying per¬ 
turbative QCD effects and in making detailed comparison with data [|10J. A variation of the 
hard scattering approach, in which transverse degrees of freedom are included in the form of 
Sudakov effects in transition form factors SB has also been employed in data analyses. 
For a critical review and comparison of the standard (Brodsky-Lepage) and modified hard 
scattering (mHSA) approaches, see Refs. |T3], [L4[]. We note that either of these approaches 


combined with data [15|] constrains the input wave-function for the quark-antiquark part of 
the //-meson. However, the gluonic part of the r]'- meson wave-function is not directly mea¬ 
sured in these experiments and will be better constrained in future experiments sensitive to 
the 7]'g*g* vertex. 

We have used the hard scattering approach to study the r] l g*g* vertex, incorporating the 
information on the wave-function and the mixing parameters entering in the g — 7 / complex 
from existing data involving the electromagnetic transitions. The results and details of the 
calculation have been presented by us in Ref. |16|. Prior to our work, Muta and Yang (TT 


derived the g'g*g vertex, with one off-shell gluon in the time-like region, in terms of the 
quark-antiquark and gluonic parts of the rf wave-function, taking into account the evolution 
equations obeyed by these partonic components. In Ref. |H|, we also addressed the same 
issues along very similar lines. We first rederived the g'g*g vertex, pointing out the agreement 


and differences between our results and the ones in Ref. |17|. The differences have to do 
with the derivation of the leading order perturbative contribution to the gluonic part of 
the 7]'g*g vertex, and the use by Muta and Yang [[17] of the anomalous dimensions derived 
in Ref. ITS] in the evolution of the wave-functions. For the anomalous dimensions, we use 


instead the results derived in Refs. 19, 201, which are at variance with the ones given in 
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Ref. ||T 8 [, but which have been recently confirmed by Belitsky and Muller |R]. Making use of 
the r]' 7*7 data to constrain the fy-meson wave-function parameters, we find that the gluonic 
contribution in the fy-meson is very significant. We then extended our analysis to the case 
when both the gluons are virtual, having either the time-like or space-like virtualities. We 
studied the effects of the transverse momentum distribution involving the constituents of the 
fy-meson and took into account soft-gluon emission from the partons by including the QCD 
Sudakov factor, following techniques which were introduced in studies of the electromagnetic 


and transition form factors of the mesons [14, 22, 231. We showed the improvements in 
the tj'g*g ^ vertex function due to the inclusion of the transverse-momentum and Sudakov 
effects, which is particularly marked in the space-like region, improving the applicability of 
the hard scattering approach to lower values of Q 2 . These effects have a bearing on the hard 


scattering approach to exclusive non-leptonic decays |24|]; the importance of the transverse- 
momentum and Sudakov effects in the decays B —> mr and B —> Kit has also been recently 
emphasized in |25|]. We also derived approximate formulae for the T)'g*g ^ vertex [16], which 
satisfy the axial-vector anomaly result for on-shell gluons and have the asymptotic behavior 
in the large-fy 2 domain, determined by perturbative QCD. 

We summarize in this report the input, salient features of the derivation of the rj'g*g^ 
vertex function, and some numerical results presented by us in Ref. [ IT> |, making comparison 
with the earlier work on this subject. 


II. DEFINITION OF THE g'g*g* VERTEX 

The ?/-meson, not being an SU(3)i? flavor-octet meson state, has a gluonic admixture 
in its wave-function in addition to the usual quark-antiquark content. We take the parton 
Fock-state decomposition of the fy-meson wave-function as follows: 


W >= sin 0 fy' > + cos 0 | 7 /' > +|rf >, 


( 1 ) 


where the SU(3fy symmetry among the light u, d and s quarks is assumed. Thus, rf q > ~ 
\uu+dd > /y/2 and \g' s > ~ |ss > are the quark-antiquark Fock states, where 0 = 39.3°±1.0° 
is the mixing angle p 6 | , and \r)' g > ~ | gg > is the two-gluon Fock state. All other higher Fock 
states are ignored. Following Refs. JI]|, |li| |20|, we express the wave function of the fy-meson 
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in terms of the eigenfunctions of the quark-antiquark |qq > and gluonic |gg > states: 


T = C[^\x,Q) + ^ 9 \x,Q) 
C = V2 f q sill 0 + f s COS 0 , 


( 2 ) 

( 3 ) 


where f q and f s are the decay constants of | rj' q > and rf s >, respectively, and their present 
estimates are: f q = (1.07 ± 0.02) f n , f s = (1.34 ± 0.06) f n , where f n ~ 131 MeV is the pion 
decay constant |26|. The eigenfucntions are then calculated for a given Q 2 by solving the 
evolution equations. As is well known, the results for both the quark-antiquark and gluonic 
components can be presented as infinite series involving Gegenbauer polynomials |L8], [19], |20f . 
Including the leading and next-to-leading terms in the expansion, the functions 0^(x, Q) 
and Q) are 


0 ^ (x, Q ) = 6xx 1 1 + 
c j^ 9 \x , Q) = xx(x — x) 


6 B^ ] 




QL S (/T ( 


b { 2 9) 

15 


165. 


(?) 


'a s (Q 


, 48 
2\\ sT 


+ 55, 


(9) 


' a s (Q- 

a s m 

as(fJ-o) 


101 1 
81 


(1 — 5xx) + 


101 1 
81 


( 4 ) 

( 5 ) 


where x and x = 1 — x are the energy fractions of two partons in the r/-meson and Q 2 > 0 is 
the energy scale parameter. It is seen that in the limit Q 2 —> 00 the quark wave-function (£|) 
turns to its asymptotic form 0 as (^) = Qxx (the same asymptotic behavior as the pion wave- 
function [pj due to its quark-antiquark content), while the gluonic wave-function (^J) vanishes 
in this limit, 0^ = 0. The coefficients of the expansion of the wave-functions (|4|) and ([5]) 
are calculated by using perturbation theory and include the effective QCD coupling a s (Q 2 ), 
which in the next-to-leading logarithmic approximation is given by []27| 

47T 


ot a (Q ) = 


201 lnlnQVA 2 
lnQ 2 /A 2 


( 6 ) 


0 o ln(g 2 /A 2 ; 

where 0o = 11 — 2n//3, 0i = 51 — 19u//3, A = Aqcd is the QCD scale parameter, and 
rif is the number of quarks with masses less than the energy scale Q. In the energy region 
m 2 < Q 2 < ml, where m c and mj are the charm and bottom quark masses, respectively, 
we have used the central value of the dimensional parameter corresponding to four active 


quark flavors, A-^ = 280 MeV |27 . 
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The invariant amplitude Ad, corresponding to the effective r]'g*g* vertex, is the sum of 
the quark-antiquark At^ and gluonic Ai^ components: 

M = M {q) +M {9) . (7) 

The quark-antiquark , and gluonic components of the g'g*g* vertex func- 


v99 


tion F, 


v a a 


entering in the invariant amplitude are defined as follows [16 


Ai (q ’ 9) = -i Ftfg* g * (q\, ql, ml) S ab £ pupa £ a *£ b *q lp q 2a , 


( 8 ) 


where £ a and e b u are the polarization vectors of the two gluons and q lp , q 2a are their four- 
momenta. 

The diagrams depicting the quark-antiquark and gluonic contents of the g'g*g* vertex are 
shown in Fig. [I] and p, respectively, and lead to the following vertex functions ||16|| : 

i 


F * ) q*A < £’ ( £’ rn v') = I dx(j)^\x,Q) 


2 Nr. 


(9) 


X 


xq\ + xq 2 — xxml + ie 


— + (x x) 


1 7T f dxtj> ,3 \x,Q) 


( 10 ) 


0 


X 


xql + xq% — (1 + xx) rri^, 


xqi + xq 2 — xxm; t , + ie 


—- (X X) 


Note that the gluonic wave-function of the rj '-meson (|j) satisfies the antisymmetry condition 
4>( 9 \x,Q) = —4>^ 9 \x,Q) JT8|, |I9|, |20fl. It implies that if the relative sign in the brackets of 
Eq. (|IU| ) were a “+” (as given in Eq. (6) in Ref. |17]], and with which we differ) the gluonic 
contribution to the q'g*g* form factor would vanish identically. 

The quark and gluonic wave-functions contain both free {B^\ B^) and constrained pa¬ 
rameters, with the latter depending on the anomalous dimensions. The free parameters can 
be fitted from the experimental data, for example, from the 7 / 7*7 transition form factor. We 
take the following restrictions on the first correction to the leading order quark-antiquark 
wave-function : \B 2 q ^\ <0.1 and \p 2 J) B^ [ <0.1 in order to keep (fM\x,Q) close to its 
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asymptotic value: 4>as(x) = 6xx, in agreement with the experimental data on the 7 / 7*7 
transition form factor [[15]. Taking into account the anomalous dimension-dependent quan¬ 


tity p^ ' 1 — —1/90 from Ref. (T^j we get 1 [L/fl < 9.0. Below, we shall present the 1 ]'g*g* 
vertex function for the maximum allowed values of the non-perturbative parameters, i.e., 


(9)1 


\B ^| = 0.1 and [-R^l = 9.0. We note that there is not much sensitivity to the variation 
of the parameter on the overall rj’g*g* vertex function . Hence, we fix this parameter 
to its maximum allowed value =0.1. However, there is considerable sensitivity to the 


2(9) 


variation of the parameter B To show this we shall take \B'f > \ = 9-0 and \B < 2 , \ = 3.0. 


(9)1 


>(9) 1 


The resulting theoretical dispersion between the two cases ([R^" = 9.0 vs. [-B^l = 3.0) 
can be seen in Fig. ^|, where we show the g'g*g vertex function for an on-shell gluon (gf = 0) 
with the other having a time-like virtuality (qf > 0) in the Brodsky-Lepage approach. The 
various leading order and next-to-leading order components and the sum contributing to 
F rt / g *g(qj,0,m^) are displayed individually. The gluonic contribution (called NLG) for the 
maximum values of the free parameters is comparable to the leading quark contribution 
(called LQ), as shown in the upper plot in Fig. [|, increasing the g'g*g vertex by almost a 
factor 2 as compared to the case when only the quark content of the r/-meson is assumed. 
This may be considered as the maximum gluonic content of the 7 /-meson allowed by current 
data. Even in the more realistic case with B^ = 0.1 and B^ = 3.0, we see from the lower 
plot in Fig. [] that the gluonic contribution is not small, and also in this case it enhances 
the value of the total g'g*g vertex function (the solid curve in Fig. ]3j) at the level of few tens 
percent. 


>(9) 1 


IV. THE 77 'g*g* VERTEX FUNCTION IN THE MHSA APPROACH 

In the Brodsky-Lepage approach to form factors, the transverse momentum dependence 
of the partons in the mesons is neglected as the hard scattering (perturbative) approach 
is applicable only when the virtualities of the external particles are much larger than the 
typical value of the parton transverse momenta kjj. Including the transverse momenta of 
the partons in the meson, the perturbative expansion of the transition form factor encounters 
large logarithms of the form ln(Q 2 /k^ i ), and it becomes mandatory to sum the multiple- 


1 This differs considerably from the values used in Q . 




gluon emissions. The formalism for the soft and collinear gluon resummation was introduced 
by Collins and Soper [^8| and by Collins, Soper and Sterman |29|]. Such gluon emissions 
give rise to powers of double logarithms in each order of perturbation theory and their 
contribution exponentiates into the Sudakov function H- The Sudakov exponents are 
known both for the quark-antiquark case, from the Drcll-Yan (DY) and the deep inelastic 
scattering (DIS) processes, and for gluons from the gluon fusion into 2q, gauge, or Higgs 
boson final states. This formalism is suitable for the description of the hadronic wave- 
functions and the hadronic form factors, such as the electromagnetic and transition form 


factors of the pion fll2[ , |22|, [23 


|, H, and has also been employed in calculating the rj'g*g* 
In the modified Hard Scattering Approach (rnHSA) [^], we take the r/-meson wave- 


vertex in Ref. 16 


function in a form similar to the pion wave-function 


23| 


T (p )(^,Q,b) = 


27tC 


0 (P) (*,Q) 


exp 


xxb 2 
4a 2 


S^\x, Q,b), 


( 11 ) 


where p = q for the quark and antiquark case, and p = g for gluons; the constant C is already 
defined in Eq. ([|), b is the separation between the //-meson constituents in the transverse 
configuration space, with b = |b| often called the impact parameter, and (fM^ and have 
the form presented in Eqs. (U) and (|5|), respectively. The transverse size parameter a can 
be determined from the average transverse momentum of the //-meson. For the numerical 
analysis, the value a -1 = 0.861 GeV is used, following from the analysis of the form factor 

|. The soft-gluon emission from the quark, antiquark and 


involving the 7r- meson [22 


gluons in the //-meson can be taken into account by including the QCD Sudakov factors 
S^\x, Q, 6 ), the details of which can be found in Ref. [|1|. 

In the space-like region of the gluon virtualities, the quark and gluonic vertex functions 
are [p~6l : 

l 


F n?l q *(Q,v,rj) = -47 ra s (Q 2 ) ° 


N r A 2 


i 


x J db\ 6 a exp 
o 


o 

xx 


dx^ q \x 1 Q) 

S^\x,Q,b)K^ + \x,b A Q A ), 


_7,2 

4a 2 A 2 A 


( 12 ) 


F l?g* g *(QiU,v) = -4? ra s (Q' 2 ) A2 


C 


x 


dx^ 9 \x,Q) J db A b A exp 

o 

\v\ Ko~\x, b\Q\) ~ {x - x)u Kq + \x, b A Q A ) 


XX 




(13) 
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Here, 

K^\x,b A Q A ) = i [. K 0 (b A Q A X + (x,u,g)) ± K 0 (. b A Q A \ + (x,u,r) ))] , (14) 

where the various dimensionless parameters are dehned as: b A = bA, Q A = Q/ A, \rj\ = 
m 2 ,/Q 2 , uj — (q 2 — cfy/Q 2 ! Ko(z) is the modified Bessel function, and 


\±(x, u, rj) — - [1 + u(x — x) ± 2 xx|t 7 |] 


(15) 


The vertex function \F r} ig* g *(qf, m 2 ,) | calculated in the space-like region is plotted in Fig. |4| 
as a function of q{ for given values of q 2 for the Brodsky-Lepage case (upper figure) and in 
the mHSA formalism (lower figure). The function q 2 F g i g * g *(qf, q 2 , m 2 ,) in the two approaches 
is shown in Fig. [5|. We note the improved perturbative behaviour of the vertex functions for 
smaller virtualities in the mHSA formalism, while for larger virtualities the two approaches 
yield very similar results. 

The transition from the space-like region of gluon virtualities to the time-like one for the 
rj'g*g* vertex function is discussed in Ref. [|T6|. The final result for the quark and gluonic 
contributions in the time-like region is: 

l 


F ifg* g *{Q,u,ri) = 4na s (Q 2 ) / dx 0 (<?) (x, Q) 


mC 


2N C A 2 


1 

x J db A b A exp 


xx 


'4a 2 A 2 A 


S^\x,Q,b)H^ + \x,b A QA), 


(16) 


F g ? g * g *(Q,uJ,r)) = 4vr a s (Q 2 ) 1 -^ J dx cj) {9) (x, Q) j db A b A 


exp 


XX 


4a 2 A 2 


b\ S^ 9 \x, Q, b) 


x 


V H [ 0 } (x, b A Q A ) + (x - x)u H^’ (x, b A Q A ) 


(+)/ 


(17) 


where 


H ( 0 ± \x,b A Q A ) = 


(18) 


H { 0 2) ( b A Q A X_(x,LJ,r /)) ± (, b A Q A \_(x,uj,r ])) 

Here, Hq 2 \z) is the second Hankel function . It is interesting to note that due to the ie 
prescription of the propagators in the hard scattering part of the g'g*g* vertex function, an 
imaginary part is generated. This can be seen in Fig. || where the upper two figures show 
the real and imaginary parts of the vertex function F,f g * g *(q 2 , q 2 , m 2 ,) as a function of q\ for 
the indicated values of q 2 , calculated in the mHSA approach. In the lower two figures, we 
compare the magnitude of the vertex function in this approach with the one in the Brodsky- 
Lepage approach. A comparison shows that there is a marked difference between the vertex 
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functions calculated in the Brodsky-Lepage and the mHSA approaches, pertaining to the 
absence of the singularity at Q 2 = m 2 , [JHJ in the latter case. 


V. INTERPOLATING EXPRESSIONS FOR THE rj'g*g* VERTEX FUNCTIONS 


For the applications in various decay and production processes it is useful to find an 
approximate expression for the vertex function which is simple and can be used over a large 
domain of the gluon virtualities. We recall that Brodsky and Lepage [Q presented an 
approximate form for the 71 — 7 transition form factor which interpolates between the PCAC 
value and the QCD prediction in the large Q 2 region. Subsequently, in Ref. |35|, this form 
was extended to the case of the 7 / — 7 transition form factor. Very much along the same 
lines, a similar expression can be written for the r]'g*g* transition form factor JTBJ : 

2v/3 UD(q 2 ,v) 


n 9 9 


(<r,^) = 47 ra s (Q- 


(19) 


V3q 2 - 87r 2 f 2 D(q 2 ,u >)’ 
where the largest energy scale parameter Q 2 = q 2 for the time-like total gluon virtuality and 
Q 2 = — q 2 for the space-like one. The function D(q 2 ,uj) is defined as follows [jHj: 


D(q 2 ,uj) — fo(uj) + 


q 


fo(u) - 


UJ 


1 - 


Q 2 l 
l — 'u 2 
2 uj 


16 B. 


(■ q) ( a s(Q 


48 

81 


In 


\<Xs(n 0 ) 
1 + UJ 


l — o; 


+ 5 B. 


g 2 {u) = 


( 9 ) f a s{Q 


101 

81 


3/o(^) - 2 

60 ; 


g 2 (uj), (20) 


where u is the asymmetry parameter having values in the interval 0 < u < 1. The asymp¬ 
totic functions for the quark-antiquark / 0 (u;) and gluon g 2 (u>) cases are derived in Ref. flT 
and obey the bounds 2/3 < / 0 (o;) < 1 and 0 < g 2 (u>) <1/6. We remark that the difference 
between the asymptotic behaviour of the vertex functions in the Brodsky-Lepage and mHSA 


approaches is small [[16 


The above expression reproduces both the anomaly value and the large Q 2 asymptotics 
of the vertex functions: 

V3 


F,yg*g 


F, » „ 
1 vg 9 


*(Q 2 ,u) 

(Q 2 ,uj 




= 4 na s (Q 2 ) 


2 \ 2UD(q 2 ,u) 


( 21 ) 

( 22 ) 


Presented in Eq. 


q 2 —>00 q - 

, the interpolating function is a smooth function in the space-like region 


of the gluon virtualities but has a pole at q 2 = 87 T 2 f 2 D{q 2 , u )/\/3 in the time-like region. 
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A similar behavior for the form factor in the time-like region was obtained by Kagan and 
Petrov [7| as the result of the evaluation of the triangle diagram. 

The mHSA approach naturally removes the unphysical singularity from the vertex func¬ 
tion in the time-like region of the gluons virtualities but the vertex function gets an imaginary 
part. In this case the real part of the approximate formula interpolates between the anomaly 
value and the large Q 2 asymptotics, while the imaginary part goes to zero as Q 2 —>• 0. An 
approximate interpolating form in the time-like region of the gluon virtualities is given in 


Ref. [16 to which we refer for further details. 


VI. SUMMARY 


We have studied the g'g*g* vertex function F ri i g * g *(q‘f,q‘f,m 2 ,) in perturbation theory for 
the most general case when both gluons are virtual. The evolution equations involving 
the eigenfunctions of the quark-antiquark and gluonic components of the 7/-meson wave 
function are solved, and the input parameters in the wave function are determined using 
data on the electromagnetic transition form factor of the 7/-meson. It is shown that within 
the allowed variation of the parameters B^ and B^ of the r/-meson wave-function, the 
gluonic contribution is not small. For extremal values of the parameters allowed by data, 
B 2 ^ = 0.1 and B^ = 9.0, the gluonic correction is found to be comparable to the quark 
contribution. But, even for a lower value B^ l} = 3.0, the gluonic contribution is present at 
the level of few tens percent. Our work corrects and extends the existing results on the g'g*g 


vertex function, reported earlier in Ref. [17[] for the case of one virtual gluon in the time-like 
region. 

We find that the Brodsky-Lepage approach leads to the appearance of a singularity in 
the region of the r/-meson mass for time-like gluon virtualities. This reflects the observation 
made earlier in the literature that the Brodsky-Lepage approach to exclusive form factors 
is valid only in the asymptotic region. In conformity with this, we find that for the total 
gluon virtuality \q 2 \ = q 2 + q 2 m 2 ,, the g'g*g* vertex function has the usual asymptotic 
behaviour: F g t g * g * ~ 1 /q 2 , anticipated for the pseudoscalar meson transition form factors. 

In the mHSA approach, where the transverse momentum dependence of the hard scatter¬ 
ing amplitude as well as the transverse momentum distribution and the soft-gluon emission 
(the Sudakov factor) in the g '~meson wave-function are taken into account, the mentioned 
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singularity in the Brodsky-Lepage-approach at q 2 = m*, disappears. Also, the validity of the 
perturbative approach is extended to smaller virtualities in the mHSA formalism, though 
asymptotically the two approaches yield very similar results for the space-like gluon virtual¬ 
ities. In the time-like region, in the mHSA formalism, we find that due to the ie prescription 
of the propagators, the rj'g*g * vertex function obtains an additional phase factor in compar¬ 
ison with the space-like expression. Thus, in applications where the r[g*g* vertex appears 
with off-shell gluons in the time-like region, this perturbative phase should be included. 

An approximate expression for the rj'g*g* vertex function is presented for the case when 
both gluons are off mass shell and have space-like virtualities. This expression interpolates 
between the anomaly value of the r]'g*g* vertex and the asymptotic QCD prediction in the 
large Q 2 region. The results summarized here have obvious applications in rare R-meson 
decays B — * r/K and B —> rfX s , and in a number of other radiative and hard processes 
involving the vertex g'g*g*. 
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qj, GeV 2 


FIG. 3: The g'g*g vertex F v ' g * g * (qf, 0, m 2 ,) as a function of q\ with B^ 1 - =0.1 and two values 
of : BIf' 1 = 9.0 (upper plot) and B^ = 3.0 (lower plot) in the Brodsky-Lepage approach. The 
dashed curves are the leading (LQ), next-to-leading quark-antiquark (NLQ), and gluonic (NLG) 
components, and the solid curve is the sum. 
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FIG. 4: The r]'g*g* vertex function | F v > g * g * (qf , g 2 , m 2 ,) | for space-like gluon virtualities in the 
Brodsky-Lepage approach (upper figure) and in the mHSA formalism (lower figure), with B^ = 0.1 
and = 3.0. The legends are as follows: q% = 0 (solid curve), g 2 = — 1 GeV 2 (long-dashed curve), 
q '2 = —5 GeV 2 (medium-dashed curve), q% = —10 GeV 2 (short-dashed curve), and q% = —25 GeV 2 
(dotted curve). 
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FIG. 5: The r/g*g* vertex function q 2 F r) i g * g *(qf,q 2 ,n r ii,) for space-like gluon virtualities in the 
Brodsky-Lepage approach (upper figure) and in the mHSA formalism (lower figure), with B^ = 0.1 
and B= 3.0. Legends are the same as in Fig. [|. 
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FIG. 6: The real and imaginary parts of the rjg*g* vertex F r) i g * g * (qf, q%, m g , ) in the mHSA formal¬ 
ism (top two figures), and its magnitude \F v t g * g *(qf,q 2 ,m^,)\ (bottom left figure) with time-like 
gluon virtualities, =0.1 and BIf'* = 3.0. The bottom right hand figure shows the form factor 
in the Brodsky-Lepage approach. The legends are as follows: q% = 0 (solid curve), r/ 2 = 1 GeV 2 
(long-dashed curve), q 2 = 5 GeV 2 (medium-dashed curve), r/ 2 = 10 GeV 2 (short-dashed curve), 
and q '2 = 25 GeV 2 (dotted curve). 






